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Abstract 

In this article we focus on the propagation of a beam of particles guided by a trans- 
versely confining potential. We consider different regimes. In the classical regime, 
we describe the beam by means of a set of hydrodynamic-like equations. We apply 
this formalism in order to investigate two practical ways for increasing the collision 
rate: by using a constriction or by tilting the guide. A high enough collision rate 
is indeed the most crucial prerequisite for reaching the quantum degenerate regime 
by means of evaporative cooling. In the quantum regime, we study the propagation 
of bosonic atoms through a constriction in two opposite regimes: the collision-less 
one and the Thomas-Fermi one. 
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1 Introduction 

So far, Bose-Einstein condensates have only been produced in a pulsed mode. 
To this respect, the achievement of a continuous beam operating in the quan- 
tum degenerate regime remains a real challenge. This would be the matter 
wave equivalent of a cw monochromatic laser and it would allow for unprece- 
dented performance in terms of focalization or collimation. Recently the MIT 
group has demonstrated the periodic replenishment of an optical trap [1]. This 
technique can therefore be used in future experiments in combination with an 
appropriate outcoupler to continuously release atoms from the trap, and hence 
to obtain a truly cw atom laser. 

Another possible way to achieve this goal has been studied theoretically in [2]. 
In this proposal, a non-degenerate, but already slow and cold beam of particles, 
is injected into a magnetic guide [3,4,5,6,7,8,9,10] where transverse evaporation 
takes place. If the elastic collision rate is large enough, efficient evaporative 
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cooling will lead to quantum degeneracy after a propagation length in the guide 
compatible with experimental constraints (i.e. a few meters). The continuous 
and pulsed loading of a slow and cold atomic beam into a magnetic guide has 
been recently demonstrated experimentally [11]. However, the collision rate is 
not sufficient so far to initiate the evaporative cooling stage. 

In this article, we start with summarizing the experimental techniques already 
developed for atom guiding. To achieve conditions favorable to evaporative 
cooling one has to maximize the collision rate. For that purpose, we describe 
quantitatively two possible manipulations of a guided thermal beam: constric- 
tion and tilting of the guide. We then study the role of interactions on a beam 
of Bose-condensed atoms propagating through a constriction. 



2 Atom waveguide 

In recent years, many experimental groups have designed both macroscopic 
and microscopic guides for transporting laser-cooled atoms. In those guides, 
the confinement is provided by means of conservative forces. So far the dipole 
and magnetic forces have been used. 

In past years, atomic transport by means of hollow optical fibers has been 
achieved [12]. In this scheme atoms are guided by the optical dipole force. 
However, the heating of atoms through intensity fluctuations and spontaneous 
emission may be a problem. 

An alternative way consists in using magnetic guides. For alkali atoms, mag- 
netic forces are strong since the magnetic moment \i is of the order of the Bohr 
magneton /ib- As the thermal energy of laser-cooled alkali atoms (microkelvin 
range) is generally much less than the ground-state hyperfine splitting (1 GHz 
corresponds to 48 mK), we can assume that the field explored in the guide 
produces a linear Zeeman shift proportional to the hyperfine g-factor gp. We 
also assume that the the spin adiabatically follows the magnetic field, yielding 
an effective potential U = gFmF(iB\B\. Depending on the atom's magnetic 
sublevel m F (positive or negative), atoms are attracted towards or repelled 
from a magnetic extremum. Actually, Maxwell's equations for magnetostat- 
ics only allow for the realization of a minimum of static magnetic fields in 
vacuum [13]. Therefore, only low-field seekers can be trapped with inhomoge- 
neous static magnetic fields. As a consequence, the atoms are not trapped in 
the energetically lowest state. This effect and the three-body recombination 
are responsible for the metastability of the confined gas. 

Quadrupole guides can be produced by means of four wires or tubes equally 
spaced on a cylinder of radius R and carrying currents +/ and —I alternatively. 
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The magnetic field is quadrupolar and is well approximated by the linear form: 
B(r±) = b'r± where r±_ = (x 2 + y 2 ) 1 ^ 2 is the distance from the z axis and 
b' cc I / R 2 the strength of the magnetic gradient. In this configuration, atoms 
with low angular momentum are not stable against spin flips. To avoid this 
loss mechanism, one superimposes a bias field B along the axis of the guide. 
Note that in this case the field strength no longer varies linearly with radius 
but is given by 

B(r ± ) = ^B 2 + b> 2 r\ ~ B + ^1. (1) 

The Taylor expansion around the center yields a quadratic variation of the 
field strength. The potential experienced by the atoms is therefore harmonic 
provided that their temperature T fulfills the condition ksT <C /J>B . A con- 
striction in the potential can be created by reducing the spacing between the 
wires of the guide and/or decreasing the longitudinal bias field. 



3 Increasing the collision rate of a classical gas 

As mentioned in the introduction, the elastic collision rate r co n in the sample 
must be high enough at the entrance of magnetic guide to start an efficient 
evaporative cooling [2]. We study in this section two ways of increasing r co n 
in the classical regime for a beam of particles propagating through the guide. 
Both methods ensure the conservation of phase space density. The first one 
consists in a constriction of the transverse confining potential, while the second 
relies on gravity to slow down the beam by tilting the guide with respect to 
the horizontal plane. 

3. 1 Constriction 
3.1.1 Single particle 

Before studying the physics of a beam in the hydrodynamic regime, it is 
worthwhile to examine the one-particle motion in a constricted guide, in 
order to gain physical insight. We therefore write the equation of motion 
along z for a particle with mass m in the power law transverse potential 
U(x,y,z) = m\(z)r]_ where X(z) denotes the z-dependent strength of the 
transverse confinement, 



Equation (2) shows that the constriction induces a transfer of longitudinal 
energy into the transverse motion. For a smooth constriction, the motion is 
adiabatic [14] and consequently the quantity i?j_A _2 ^ 2+7 ^ is constant, where 
E± denotes the transverse energy. As the total energy E = E± + mz 2 /2 
is also conserved, the decrease in longitudinal velocity can easily be related 
to the increase in A. We notice that the effect of the constriction is more 
dramatic for particles having a high transverse energy, and is absent for a 
particle propagating along the axis. 



3.1.2 Hydrodynamic regime 

We now assume that there are enough collisions to ensure local thermody- 
namical equilibrium. From Appendix A we deduce the equations governing 
the beam temperature T(z,t), linear density n(z,t) and mean longitudinal 
velocity u(z, t): 
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We now study the stationary regime. The previous equations become 
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where <fi and u are constants. Eq. (3) corresponds to the conservation of the 
number of particles, Eq. (4) of momentum along z and Eq. (5) of enthalpy. 
After some algebra one obtains the relation between the strength of the poten- 
tial A = A(z)/A(0) normalized to its initial value and the reduced temperature 
f = T(z)/T(0): 



A = T^- 2 )/ 4 fl--(T -1)) 
V Vo ) 



7 /4 



(6) 



where we have introduced the parameters 5 = 5+4/7 and r] = m-u(O) 2 / (k B T(0)) 
which compares the initial kinetic and thermal energies. For typical experi- 
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mental values [11], one has f) > 1 and (6) yields 

f ~ ^4/(7(5-2)) _ ^ 

As the normalized collision rate f = r co u(-z)/r co u(0) increases as f ~ T 2 , 
one deduces the scaling law T ~ A 8/(7(5 ~ 2)) , which shows that moderate con- 
strictions (for example by a factor of 4) can increase the collision rate by 
a large amount (nearly one order of magnitude for an harmonic transverse 
confinement). 



3.2 Tilted guide 



We now consider that the guide axis z is tilted by an angle a with respect 
to the horizontal plane. The strength A of the transverse confining potential 
is assumed constant. By tilting the guide, the mean velocity of the beam 
decreases while the linear density and the temperature increase, resulting in 
an increase of the collision rate. Following the same procedure as the one 
developed in Appendix A, we get, in the stationary regime, the following 
equations : 
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where g is the acceleration of gravity. This set of equations can be solved 
analytically, one deduces the relation between the reduced temperature T and 
the reduced velocity of the beam u = u(z)/u(0): T = u 2 ^ 2 ~ s \ We also establish 
the relation between the collision rate and the coordinate z = z/z with z = 
11%/ (g sin a), and the link between the parameter f) = mu 2 /(riotcsT) = u 2 /T 
and z: 



z(f) = 6 + r]0 _ J_fV2 _ If (2-«5)/2 
2r] 2 Vo 2 

m ~ 2Vo 2^ 2 V 

We recall that for a particle with initial velocity Uq moving in this potential, 
the turning point corresponds to z = 1/2. Fig. 1 (a) shows the gain in the 
collision rate as a function of z for 2 values of the parameter i] = 100 and 
r] = oo and for an harmonic (7 = 2) and a linear (7 = 1) confinement. Along 
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the propagation through the tilted guide, the temperature T increases and 
the mean velocity u decreases resulting in a strong reduction of the parameter 
r)(z) as a function of z (see Fig. 1 (b)). 

In this section, we have described the propagation of a beam of classical parti- 
cles through a constriction and a tilted guide by means of a set of hydrodynamic- 
like equations. We have used this formalism to investigate two practical ways 
for increasing the collision rate. Indeed those results are valid for a beam 
of bosonic particles. We recall that a high collision rate is the most crucial 
requirement for reaching the quantum degenerate regime by means of evap- 
orative cooling [2]. The remainder of the paper is devoted to the quantum 
degenerate regime. We focus on the propagation through a constriction, and 
analyze the role of interactions. 



4 Constriction of a quantum beam 

For atomic waveguide geometries with very small transverse dimensions as now 
experimentally demonstrated [3,4,5,6,7,8,9,10], the quantum nature of atoms 
starts to influence the dynamics. This problem has first been investigated 
in the context of quantized-conductance in a two-dimensional electron gas 
[15]. In those models, the conducting properties of a nanoscopic constriction 
are investigated in the free-electron approximation with hard-walls boundary 
conditions. In this section, we address the corresponding problem for a beam 
of Bose-condensed atoms confined in the transverse direction. In the following, 
we restrict our study only to harmonic confinement for the sake of simplicity. 
As for the classical case of the previous section, we model the constriction by 
a dependence of the transverse frequency with the longitudinal coordinate. 
The potential experienced by the atoms is therefore of the form U (x, y, z) = 
muj(z) 2 [x 2 + y 2 }/2. We begin with the case of an ideal beam, and then consider 
the role played by interactions for the propagation through a constriction. 



4-1 The ideal Bose-condensed beam 

In this paragraph we study the propagation of an ideal Bose-condensed beam. 
For this purpose, we solve the stationary Schrodinger equation by means of 
the adiabatic basis (the harmonic oscillator's one with an angular frequency 
that depends on the longitudinal coordinate z). Details of this lengthy but 
straightforward calculation can be found in the Appendix B. If k denotes the 
incoming longitudinal wave vector, the local wave vector fc ,o for atoms in the 
transverse ground state within the adiabatic approximation reads (see Eq. 
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(B.2)): 



kr. 



0,0 — 




(11) 



a being the local harmonic oscillator length <r(z) = (h/muo(z)) 1 / 2 . The physical 
picture is quite clear: the constriction induces an augmentation in the zero- 
point energy for the transverse degrees of freedom. As the total energy of 
the atom is conserved, the longitudinal kinetic energy decreases. Actually, 
one recovers, in the quantum domain, the transfer of kinetic energy from 
longitudinal to transverse degrees of freedom as already pointed out for the 
classical dynamics. However, the constriction modifies the propagation even 
for particles in the transverse ground state because of the finite size of the 
transverse wave function. 

4-2 The role of interactions 

We investigate in this paragraph the role of interactions on the propagation 
through a constriction of a beam of interacting Bose-condensed atoms. The 
main difficulty arises from the non-linearity of the Schrodinger equation due to 
the mean field interaction term. The strength of the interaction is proportional 
to the scattering length a [16]. By assuming that the transverse dimension of 
the cloud has the local equilibrium shape adapted to the local number of 
particle per unit length, the problem becomes one-dimensional and one can 
introduce a solution for the wave function of the form ^(r, t) = f(z, t)g(x, y; n) 
in the variational procedure [17], where n is the local density of particles per 
unit length. By writing the longitudinal part of the wave function / in the 
phase and modulus form / = y/ne^, we introduce the longitudinal velocity 
field defined by u = (h/m)dtp/dz. 

After some algebra (see Appendix C), one gets the stationary hydrodynamic 
equations for the propagation of a beam of interacting Bose-condensed parti- 
cles: 

nu = (f) (12) 



where ft, is an effective chemical potential, and and c two constants depending 
on the initial conditions. Eq. (12) reflects the conservation of the flux while 
Eq. (13) accounts for the conservation of energy. 

As a first example, from Eq. (12) and (13) we recover the results of the previous 
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section for a perfect gas for which a = and jl = fiio in a smooth constriction 
i.e. such that the last term of the l.h.s of (13) is negligible. In the following, 
we apply this set of equations in the weak-coupling limit (na <C 1) and in 
the Thomas- Fermi limit (na 3> 1). For those two limits, the value of (1 differs 
(see Appendix C). It is convenient to introduce the dimensionless notations: 
u = u(z)/u(0), z = z/a(0), Co = u(z)/u(0), a = mw 2 (0)/(fiw(0)) and f3 = 



y/na = \/a<fi/u(0). The combination of Eq. (12) and Eq. (13) leads to the 
relation between the reduced velocity u and the rate of compression Co as a 
function of the dimensionless parameters. 

In the weak-coupling limit that corresponds to the low density regime (/? <gC 1) 
we find: 

d 2 ii 3 f dii\ 2 

^ = — 1 — 1 -2au(u 2 -l)+4u(l-u) + 8u[3 2 (l-u;/u) (14) 



while in the opposite regime, commonly called the Thomas-Fermi regime, we 
have 



In the following, we consider a linear increase of the transverse confining fre- 
quency: Co = 1 + 72. The results of the numerical integration of Eq. (14) and 
(15) are plotted on Fig. 2, for the following values of the dimensionless pa- 
rameters: a = 10 4 , 7 = 10" 4 , (3 having the value 10 in the Thomas-Fermi 
regime, and 0.1 in the weak-coupling regime. The large constriction ratios im- 
plied by 7 = 10~ 4 can be achieved on microchips [18]. For a given compression 
rate, the reduction of the velocity for the propagation through a constriction 
is dramatically magnified in the interacting case. It reflects the fact that the 
ground state wave function has a much broader transverse profile because of 
repulsive interactions (y/na 3> 1). As already emphasized in §3.1.1, the larger 
the transverse size of the beam, the higher the coupling between longitudinal 
and transverse degrees of freedom. We point out that this magnification of the 
coupling between degrees of freedom induced by interactions may have to be 
taken into account for interferometers based on guided matter waves [19]. 



5 Conclusion 



In this article, we have investigated quantitatively two possible ways (by com- 
pression and tilting of the guide) of increasing the collision rate of a thermal 
guided beam while keeping the phase space density constant . Those results 
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are important from an experimental point of view since as soon as the collision 
rate is high enough, one can implement evaporative cooling as suggested in [2] 
to reach the quantum degeneracy regime. Those results have been obtained 
under the hydrodynamic formalism. In future work, we intend to the study the 
case where the collision rate is not sufficient to ensure the local equilibrium 
requirement for hydrodynamics. We have also carried out the calculation for 
the propagation of a beam of bosons through a constriction in the degenerate 
regime in presence or not of interactions between particles. Strictly speaking, 
those calculations are valid at T = 0. Actually, the propagation of a beam of 
bosons through a constriction below the critical temperature can be regarded 
as a two-fluid problem. To address this problem completely, one would have 
to take into account the interaction between the thermal and the condensed 
part. 
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A Theoretical framework for guided beams in classical regime 

We consider a gas of particles with mass m confined by a static, velocity- 
independent potential U(r). In this appendix we derive the classical equations 
for a beam of particles confined in a cylindrical geometry. For that purpose, 
we start with the classical Boltzmann equation: 



where /(r, v, i) is the distribution function and 7 C oii[/] is collision term. The 
density p is obtained from the distribution by integration over velocity: 



dt dr m dr dv 



(A.l) 




(A.2) 



9 




Fig. 1. (a) Collision rate normalized at its initial value f (z) = r co n(£)/T co n(0) as a 
function of the reduced longitudinal coordinate z for a tilted guide, for rjo = 100 and 
7/0 = oo. The dashed curves correspond to a linear transverse confinement (7 = 1) 
and the solid lines to a harmonic transverse confinement (7 = 2). (b) Parameter 
7/ = mu 2 (z) I (ksT(z)) as a function of z, for an initial value of 100, for linear (dashed 
line) and harmonic (solid line) transverse confinement. 
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Fig. 2. Longitudinal velocity u = u/uq as a function of the constriction ratio Co 
for a quantum beam in the weak-coupling regime (dashed line, (3 = 0.1) and in 
the Thomas- Fermi regime (solid line, (3 = 10) for a = 10 4 and 7 = 10 -4 . In the 
Thomas-Fermi regime, the repulsive interactions lead to a strong coupling between 
longitudinal and transverse degrees of freedom resulting in a larger reduction of the 
beam velocity through the constriction. 

A.l Reduction of dimensionality 

We present the formulation of conservation equations for cylindrical geome- 
tries. We choose the longitudinal axis along z. In order to get hydrodynamic- 
like equations, we assume that the distribution function / is given by a local 
Maxwell-Boltzmann distribution of the form: 



In the ansatz (A.3) we assume that the velocity field is purely longitudinal, 
and that all macroscopic quantities depend only on z [20]. We introduce the 
linear density n defined by: 



The conservation of the number of particles is obtained by taking the av- 
erage of (A.l) over the transverse coordinates and the velocity components: 




(A.3) 




(A.4) 
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/[A.1] dxdyd 3 v. 




(A.5) 



where (• • •) denotes the following normalized average 

(A) = ^- [ dxdyd 3 vAf(r,v,t). 
n{z) J 

The axial momentum current equation is deduced from the average of the 
on-axis component of the velocity J[A.l]i> z dxdyd 3 v: 




(A.6) 



The conservation of energy is obtained from the average of the one-particle 
hamiltonian H(r,v) = mv 2 /2 + U(r), J[A.l]H(r, v) dxdyd 3 v 

Fti n{H) ) + iL{ n{vz){H) ) + 7 £( n(( ^ 3) " K)( ^ 2)) ) = (A - 7) 



A. 2 Application to a power law potential for a stationary flow 



We consider in this section that the transverse confining potential is given by 
a power law of exponent 7 with a z-dependent strength \{z): 

U(x,y,z) =m\(z)rl (A.8) 



where r±_ = (x 2 + y 2 ) 1 ^ 2 denotes the radial distance from the guide axis. 

Furthermore, we define the mean velocity of the beam by u = (v z ), and calcu- 
late with the ansatz (A. 3) the following useful quantities: (v 2 ) = u 2 + ksT/m 
and (v 3 ) = u 3 + 3uksT/m. One also readily calculates the average values 
(U) = 2k B T/-f, (H) = (2/-f + 3/2)k B T + mu 2 /2. The one-dimensional hydro- 
dynamic equations (A.5) to (A. 7) then simplify considerably. In the stationary 
regime that we shall consider hereafter, they can be recast in the form: 
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An interesting feature of this set of equations is the conservation of phase-space 
density D = pA|g, where X^b — h/ \f2nmkBT is the thermal de Broglie wave- 
length and p the 3D particle density. Indeed, for the power-law potential (A. 8), 
D scales with the four parameters u, T, A and 7 as: D ~ A 2 / 7 -u _1 T _2 / 7_3//2 . 
One can easily check that the set of equations (A. 9) implies that 

dD 

7i7 = - 



The collision rate r coll scales as pT 1 / 2 , thus the ratio f = r co n(z)/r coll (0) scales 
as (T(,)/T(0)) 2 . 



B Compression of an ideal quantum gas 



In this appendix, we address the problem of the propagation of a quan- 
tum gas in a harmonic waveguide with a constriction. The potential expe- 
rienced by the atoms reads: U(x,y,z) = muj 2 {z)[x 2 + y 2 )/2, cu(z) being the 
z-dependent angular frequency. We introduce the harmonic oscillator length 
a(z) = [h/mujiyz)) 1 ^ 2 , and denote its first and second derivatives with respect 
to z by a' and a". For a theoretical description, we resort to the study of the 
stationary solutions of the Schrodinger equation [15]. A general wave function 
can be expanded on the adiabatic basis as \^/) = J2n x ,n y a n x ,n y {z)\n x (z)) ® 
\n y (z)), where (|nj(,2))} denotes the local ID harmonic basis for the i-direction 
with a z-dependence in angular frequency. Inserting into the Schrodinger 
equation with stationary energy E = h 2 k 2 /2m and projecting on the mode 
\m x (z)) <S> I m y (z)) yields the following set of coupled linear equations: 



^ z 2 ~ Xft"™!* y ^ z ^rn x ,n x > ,m y 'm, ,n x J 



d ( \ ni , .111 \ / .. * l-O , 

Ef da mxt n y 
\ z -J- i \,n, ' a m x ,ny-Dm y ,n y 
n y ^m y V UZ 

+2 oz nxtriy A mxtnx A myiny , (B.l) 

where we have defined the matrix elements A muni = {mi\d/dz\rii), B muni = 
(rrii\d 2 /dz 2 \ni), and the local wave number k mx)Tny 



a' 2 . 9 9 . 1 + 2m x + 2m v 
k 2 mx , my = k 2 -—(2 + m x + m y + m 2 x + m 2 y ) y -. (B.2) 

ZiCJ u 
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The matrix elements can be readily calculated: 

A mi , rH =a'{m,i\af - a 2 i \n i )/2a 

B mi , ni = -^(o' 2 {mi\[af - a 2 } 2 \m) - 2(aa" - a' 2 ){mi\af - a 2 [n^, 

where a\ and Oj are the creation and annihilation operators for the considered 
harmonic basis. 

Eq. (B.l) shows that each mode experiences a different potential barrier. In 
particular, k mx . my (z) decreases in for a propagation through a constriction 
since the ground state energy increases in this case. From the above analy- 
sis we can conclude that an incident wave initially in the transverse ground 
state remains in this state if the adiabatic approximation criterium is fulfilled, 
namely \a'\ <C 1. The adiabatic approximation fails around the turning point 
if the constriction is such that the incident wave can be reflected. We note that 
if the condition of adiabatic approximation is not fulfilled, Eq. (B.l) allows 
one to evaluate the importance of inter-level scattering. 



C Effective one-dimensional hydrodynamic equation in the quan- 
tum regime 

In this appendix, we derive the effective ID hydrodynamic equations for a 
beam of Bose-condensed interacting atoms propagating through a guide. For 
this purpose, we recall the main steps of their derivations as presented in [17] 
by means of the standard variational procedure. This approach is valid under 
the dilute-gas assumption [16], namely pa 3 <C 1 where a is the scattering 
length and p the mean density. 

C.l Hydrodynamic equations 

The action S for the condensate wave function ip reads: 

where the energy functional £ is given by 

£[W\= J d 3 r(^-|V^| 2 + ^|^| 4 + f/(r)|^| 2 ). (C.2) 
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The fact that the transverse beam size is small compared with the other rel- 
evant length scales allows us to consider the density profile across the guide 
to have its equilibrium form appropriate to the local number of particles per 
unit length. As a consequence, the problem becomes essentially one dimen- 
sional. In this case, the wave function can be assumed to have the form: 
i/)(r,t) = f(z,t)g(x,y;n) where n is the local density of particles per unit 
length and g the equilibrium wave function for the transverse motion. We 
emphasize that this ansatz does not correspond to a complete separability 
between the axial and transverse coordinates since g depends on z via the 
density per unit length n. In the following, we choose g to be normalized to 
/ \g\ 2 dxdy = 1, and therefore \f\ 2 = n. 

The equations for g and / are obtained by minimization: 

^(£-^/d 3 rH 2 )=0 and ^ = 0, (C.3) 



where /x is the chemical potential. After some algebra, one obtains the following 
set of equations: 



-— A + U+ \fg\ 2 )9 = H9 
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Anh 2 a 



m 



f\f\ 2 J dxdy\g\ 4 + fj dxdyU\g\ 



(C.4) 



(C.5) 



where we introduce the effective chemical potential p, = fi—h 2 \df/dz\ 2 / (2m|/| 2 ) 
To obtain the hydrodynamic-like equations, it is convenient to write / in terms 
of a modulus and a phase: / = y/ne' 1 ^. The velocity field u associated with 
the phase f reads u = (h/m)dip/dz. From (C.4) and (C.5), one derives the 
equations for the density n and the velocity u: 



du d /_ i 2 h 2 i 9 2 (v^)\ 

m m=-dz-{ fI+ 2 mu -*^rd*-y (c - 7) 

The first equation is the usual equation of continuity, while the second is an 
Euler-like equation. In the adiabatic limit, the last term of (C.7) is negligible. 
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C.2 Determination of the chemical potential 

In the following, we derive the expression for the chemical potential in two 
opposite limits, namely the Thomas Fermi regime (na ^> 1) and the weak 
coupling limit (na <C 1). We restrict our analysis to an harmonic trapping 
potential: U(x,y,z) = mu 2 (x 2 + y 2 )/2. 

The Thomas-Fermi regime 

In this large density limit (na ^> 1), the kinetic term in (C.4) can be neglected. 
The normalization condition on g provides the relation between the effective 
chemical potential ft and the local density of particles per unit length: ft = 
2fko\/na. If u) does not depend on z or changes smoothly, \df/dz\ ~ and 
ft — //, we therefore recover the results of Refs. [17,21]. 

The weak-interaction limit 

In the absence of interaction, go has a Gaussian form, |go| 2 = (no- 2 )~ l e~ <yx2+v2 ^ a ' 



with a = yh/mu and /i = h~u. The low density limit corresponds to na 1, 
we consequently derive perturbatively from Eq. (C.4) the correction to the 
chemical potential: 
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